Fluctuation theorems for an asymmetric rotor in a granular gas 
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We investigate the validity of fluctuation theorems for an asymmetric rotor experiment in a 
granular gas. A first state, with a Gaussian distribution of the angular velocity, is found to be 
well described by a first order Langevin equation. We show that fluctuation theorems are valid 
for the injected work and for the total entropy production. In a second state the angular velocity 
distribution is double-peaked due to a spontaneous symmetry breaking: A convection roll develops 
in the granular gas, which strongly couples to the rotor. Surprisingly, in this case similar symmetry 
relations hold, which lead to a good prediction for the height ratio of the two peaks. 
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Under rather general conditions, the thermodynamic 
description of non-equilibrium systems is constrained by 
the consequences of the fluctuation relations In par- 
ticular, the fluctuation theorem (FT) provides, for a sys- 
tem in a non-equilibrium steady state, a quantitative 
symmetry relation between the probability of having a 
positive fluctuation for the entropy production in a time 
t and a corresponding negative one: 

^( ^ T=+<l l ) =P a Va forr»r c (1) 
\P(a T = -a)J 

where r c represents the largest characteristic time of the 
system and f3 a prefactor [2-0, EH . In general, the hy- 
potheses used to prove this theorem are not verified ex- 
perimentally and therefore, it is not clear whether Eq. (fTJ) 
holds or not. Experimental tests of FTs have been mostly 
performed on stochastic systems in contact with a ther- 
mal bath. In such experiments, a FT is valid for the 
injected work into the system or the total entropy pro- 
duction and the coefficient f3 is directly related to the 
temperature of the thermal bath j3 = 1/fcsT 

Experimental tests are particularly scarce and incon- 
clusive for another type of systems for which FTs are be- 
lieved to hold, namely non-thermal dynamical systems. 
In these systems, the strong fluctuations come from the 
non-linear interaction of many degrees of freedom of a 
dissipative system. An interesting example of such a sys- 
tem is a granular gas, which by its nature is already out- 
of-equilibrium due to the dissipative character of the in- 
elastic collisions. Experiments searching for FTs in gran- 
ular systems have been performed [|| , but their interpre- 
tation remains unclear [5j . It is a challenge to study FTs 
in such a system and to test to what extent the fluctua- 
tion relations can be satisfied. 

In this letter we will study FTs for a rotor immersed 
in a granular gas. Eshuis et al. 0] showed that such a 
system not only exhibits Brownian-like dynamics com- 
parable to a thermal system, but also a state in which 
symmetry is spontaneously broken, and for which there 
exists no thermal analogue. Here, we will show that FTs 
are relevant to both states. 



Experiment - The rotor is composed of four vanes 
(25 x 60 mm 2 each, made from a single piece of stain- 
less steel), precisely balanced around an axis which is 
connected to the container by a low-friction ball bear- 
ing. The granular heat bath consists of glass beads of 
diameter d = 4 mm (density p = 2600 kg/m 3 ), which 
are fluidized by vertical vibrations of the bottom (z = 
at rest) with amplitude a and frequency / such that the 
grains interact dissipatively with the vanes. We present 
the results obtained for a system of TV = 1000 parti- 
cles (500 and 2000 particles have also been used) and 
for two different heights of the axis above the bottom 
(h = 51 mm and h = 75 mm). The natural dimension- 
less control parameter of the granular heat bath is the 
shaking strength S = Air 2 f 2 a 2 /(gh), which represents 
the ratio of the typical kinetic energy injected into the 
system by the vibrating bottom and the potential energy 
of the particles at the height of the axis of the rotor. 
The granular temperature T g is denned as the mean of 
the kinetic energy fluctuations per particle. The angu- 
lar position of the vanes is measured using an optical 
angle encoder, with an accuracy of 1.9 • 10~ 7 rad, at a 
sampling rate of 1000 Hz (larger than the typical colli- 
sion rate, which is about 100 Hz). After the system has 
reached a steady state, we start to record the position for 
typically 15 min, which is about 10 3 times the relaxation 
time of the system. We repeat experiments several times 
with the same conditions (S and h). An asymmetry is 
voluntarily introduced in the system by coating the left 
hand side of each vane with rubber tape. This side will 
therefore be softer, diminishing its coefficient of normal 
restitution. As a results, vanes are expected to preferen- 
tially rotate in counter-clockwise direction (positive 0). 
Such a system behaves like a ratchet, as explained in |7J- 

Q. 

The dynamics of the system strongly depends on the 
experimental parameters (height of the vanes h, shaking 
strength S, and number of particles). Fig. [T^, presents 
two time series of the angular velocity u = 6 of the rotor: 
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State I has an angular velocity which fluctuates around 
a non-zero positive mean value (due to the symmetry- 
breaking coating). In state II, the angular velocity is 
either positive or negative and there are switches from 
one state to the other one. This is a state in which, in 
addition to the symmetry-breaking coating, there is a 
spontaneously broken symmetry, which can be explained 
from the development of a convection roll in the granular 
gas stabilizing the motion of the vanes 0. These two 
states will now be studied in detail. 

State I - In the first state, the dynamics of the angular 
velocity uj can be decomposed into a constant positive 
average value {uj) and strong fluctuations Suj = uj — (uj) 
around this mean. To shed more light on the dynamics 
and statistics, the power spectral density (psd) of Suj and 
the probability density function (pdf) of uj are plotted in 
Figs. [Tp and c. The Lorentzian shape of their psd and 
their Gaussian distribution indicate that the fluctuations 
Su> are identical to those in the absence of the ratchet 
effect, as expected for a linear situation. The dynamics 
is therefore likely to be well described by a first order 
Langevin equation 



r du) 



-juj + M 9e + i] , 



(2) 



dt 

where / (= 7.2 • 10~ 5 kgm 2 ) is the angular moment of 
inertia, 7 a viscous drag coefficient, and r\ a stochastic 
noise, ^-correlated in time, due to the random collisions 
between the particles and the rotor. This description 
is compatible with kinetic theory in the limit in which 
the mass of the ratchet is much larger than the mass of 
the particles 
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The measurement of the variance of 
the fluctuations, T T /I, and the cut-off frequency of the 
psd, / c = 7/(271-/) (related to the relaxation time of the 
system, to = l/(27r/ c )) provides the value of the viscos- 
ity 7 w 7.0 • 10~ 5 kgm 2 /s and the ratchet temperature 
T r « 0.115 mJ, which is close to the granular tempera- 
ture [l7| ■ The ratchet effect is described by adding a con- 
stant torque Mg e to the Langevin equation llj , which is 
identified with 7(0;) (« 4.5 ■ 10 -5 Nm) by time averag- 
ing of Eq. (T2J) . The subscript de points to the fact that 
|Ma e | is an increasing function of |<9e|, the difference in 
normal restitution coefficient of the coated and uncoated 
side of each vane. For the symmetric system, i.e., without 
coating, (ui) = and thus Mg e = 0. 

Based on the Langevin description, the amount of 
work injected into the rotor during a time r is W T = 
J t t+T M de uj(t')dt' = M de A6, where A6 = 9{t + r) - 0(t), 
i.e., the symbol A from hereon represents the difference 
measured over a time delay r. For W T we expect a FT 
of the form 

P(W T ) \ . ( P(A0) 



In 



P(-W T ) 



In 



P(-A9) 



M de A6 



(3) 

where P{A0) is the pdf of AO. The expression in the left 
hand side of the equation is called the symmetry func- 
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FIG. 1. (color online) (a) Typical time series of the angular 
velocity uj for state I [S = 2.78 and h — 75 mm (top, red)] 
and for state II [S = 2.15 and h = 51 mm (bottom, blue)], 
(b) Corresponding power spectral densities of the fluctuations 
Suj, which for both cases virtually overlap and are fitted well 
by a Lorentzian. (c) Corresponding probability distribution 
functions of uj. For state I (red symbols) the pdf is well fitted 
by a Gaussian (red line). For state II (blue symbols) there is a 
pronounced double peak; the blue line is the result of Eqs. (j6j) 
and (Q. 



tion. Due to the Gaussianity of the distribution of w, the 
pdfs of AO arc gaussian too, and the symmetry functions 
are expected to be linear with AO for different values of 
t. This is indeed found in experiment (Fig. [5Ji). The 
slope, £(t), is found to decrease with r and reaches the 
constant value Mg e /T Y expected from the FT Eq. ([3]) for 
t > 3 s, which is approximatively 3 times the relaxation 
time of the system. Thus we find that two independent 
measurements -one from the pdf of Fig.QJi and one using 
the FT Eq. ([3])- lead to the very same value of the rotor 
temperature T r , hereby confirming the validity of the FT 
Eq. ©. 

We now turn to the fluctuations of the total (trajec- 
tory dependent) entropy production in a time span r, 
Astot, as defined by Seifert in [ljj. We follow Q and 
first write the dissipated heat Q T [l8| as the injected 
work W T minus the difference of the kinetic energy of 
the rotor Ek = \Iuj 2 at the beginning and end point. 
The entropy change As m in a time span r is now defined 
as the dissipated heat divided by the ratchet temperature 



T r As m = Q T = W T - AE, 



k ■ 



(4) 
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FIG. 2. (color online) State I [S = 2.78 and ft = 75 mm]: (a) The symmetry function ln(P(A0)/P(-A0)) versus A0 for different 
values of the time interval r. All are linear with A9, such that we can compute their slopes S(r). (b) S(r) as a function of r. 
In the limit when r is larger than the relaxation time of the system E(r) is equal to Mg e /T r (dashed line), as expected from 
the steady state FT Eq. Q. (c) The symmetry function ln(P(As to t)/-P(— As to t)) versus the total entropy production As to t 
for different values of r. As predicted by the detailed FT, the data lie on a straight line with slope one (dashed line). In (a) 
and (c), the deviations from the linear regime for large fluctuations are due to a lack of statistics in a region where the pdfs 
are at least two orders of magnitude smaller than their maximum value. 



First we note that for the symmetric system (i.e., without 
the symmetry-breaking coating) the dissipated heat is 



equal to AE, 



k,sym 



= T r Ai 



m,sym 7 



which is small but not 



strictly zero. This contribution has to be subtracted from 
the entropy change Eq. ^ in order to obtain the entropy 
created by the presence of the external torque Mg e . As 
argued in Q , it is this last entropy that may be identified 
as the total entropy production, i.e., As to t = As m — 
As mjSym . This quantity is interesting because it should 
satisfy a detailed FT: The symmetry functions for Astot 
are expected to be equal to As to t itself for all values of 
r when the system is in a steady state lU 14 1. From 
the Langevin equation l[2]). the total entropy takes the 
following form 



T r As to t = M de A0 - /(oj)Aw. 



(5) 



In Fig. [5t we present the symmetry functions for dif- 
ferent integration times r: We find that the pdfs are 
Gaussian (not shown) and, most importantly, that the 
symmetry functions lie on a straight line with slope one 
within experimental errors, whatever the time delay r 
is. So indeed a detailed FT holds for the total entropy 
production. 

State II - In state II the dynamics is completely dif- 
ferent due to the coupling to a spontaneous convection 
roll in the granular heat bath: The time series of Fig. [TJi, 
shows fluctuations around two preferred velocities, cor- 
responding to the two rolls. As a result, in the pdf of 
w (Fig. [Tfc) two distinct peaks appear. Their different 
heights reflect that the system obtains a finite mean an- 
gular velocity (ui) by residing preferably around the pos- 
itive peak - in contrast to the shift of the entire pdf in 
state I. From Fig. [TJd we find that in both states the psds 
are almost exactly the same and very well fitted by a 
Lorcntzian. This indicates that the fluctuations in each 
of the two peaks -and therefore the short time dynamics- 



are very similar to those of state I. 

This suggests the separation of the dynamics into two 
parts: On a short time scale, oj fluctuates around a mean 
value, with dynamics described by a Langevin equation 
similar to Eq. ([2]), independent of the direction of the 
roll. On a much larger time scale, a coupling between the 
granular gas and the vanes induces collective motion, in 
which every now and then, through a particularly strong 
fluctuation, the mean value of uj quickly switches between 
positive and negative, reversing both the sense of rota- 
tion of the vanes and the roll. These reversals happen 
randomly in time. The system stays few seconds in each 
state as can be seen in Fig. QJ,. Due to the asymmetric 
coating of the rotor, reversals are easier to realize when 
the vanes are rotating in the clockwise direction than 
anti-clockwise, explaining why the system has a prefer- 
ence for the latter 

Guided by the above considerations, we write the pdf 
of u> as the sum of two Gaussian distributions with the 
same variance a 2 = T r /I; one with mean value u>o + ujg e 
and the other with — ujq + ujg e : 



P(u) = A 



P+ exp ( - — (w - (ujq + uJde) 2 



(6) 



Here ujQ e stands for the shift due to the symmetry 
breaking coating (i.e., ±wo denote the locations of the 
peaks without coating). The factors p± represent the 
weight of each peak (the probability of each direction of 
the convection roll) and the normalization factor A = 
I j '{iirT-c) 1 / 2 has been chosen such that p+ + p~ = 1 im- 
plies J P(uj)doj = 1. 

What can we learn from FTs for this very non- 
Gaussian system? It seems reasonable to assume that if a 
FT exists, it will be for the work done on short timescales, 
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FIG. 3. (color online) State II [S = 2.15 and ft = 51 mm]: 
(a) Symmetry function of A8 for three different integration 
times, (b) Symmetry function of u>. 



when the system resides in either of the two peaks of the 
pdf of u>. And indeed, in Fig. [3^ (dashed blue line) we 
find that the symmetry function for AO is linear in AO 
for t = 0.01 s. Since in this limit uj ~ A(9/r, we may 
expect a similar symmetry relation for uj, which we plot 
in Fig. \3]p. Indeed, In \P(uj) / P(— uj)] oc uj, which is truly 
remarkable in view of the peculiar shape of its distribu- 
tion (Fig. QJ). Computing the symmetry function from 
the double Gaussian of Eq. (j6]) we find that it can only 
be linear if the weights of the peaks obey 



P+ (2IuJ de UJ 
P- V T r 



which then directly leads to the symmetry relation 



In 



P(-uj) 



2Iujg e 



-UJ . 



(7) 



(8) 



From the slope of the experimental symmetry function 
(Fig. [5}d), which, with T r = 0.062 mJ leads to ujq & = 
0.372rad/ s rad/s, combined with ujq = 1.8 rad/s from 
the pdf of the symmetric system, we compute p+ and 
p_ and plot the resulting theoretical pdf in Fig. QJ. The 
agreement with the experimental pdf is good, except for 
angular velocities close to zero. 

We now return to the symmetry function of A9 in 
Fig. [3ja,): For the smallest value (r = 0.01 s) we find 
a symmetry relation. This is because the changes AO are 
dominated by the fluctuations, and the reversals are just 
rare events without large consequences for the value of 
A8. For the intermediate value r = 3 s, there is no valid 



symmetry relation. This stands to reason because when 
the integration time r increases, the probability that it 
contains at least one reversal increases as well. AO is 
then likely to include a reversal and there is no reason to 
consider the system to be in a steady-state, which is a 
prerequisite for a symmetry relation to hold. There are 
two distinct regimes: One at small \A0\, dominated by 
reversals, and another at large \A0\ > uj t, a distance 
which in the given time interval is unlikely to be reached 
when the trajectory includes one or more reversals (dot- 
ted red line in Fig. [3^,). For very large r, typically many 
reversals are included in each AO and we can consider 
the system to be in a steady state again. Indeed, for 
t = 10 s the pdf is Gaussian (not shown) due to the cen- 
tral limit theorem, and the symmetry function tends to 
become linear (Fig. solid black line). However, the 
experimental slope is smaller than the slope Mg e /T r one 
would expected based upon Eq. ([3]). This may be be- 
cause the entropy production depends on two parts: one 
is the work Mg e A0 and a second one coming from an- 
other variable responsible for the reversals, which is not 
included here (cf . the theoretical work in [15j ) . 

In conclusion, we have investigated the validity of FTs 
in two different states of an asymmetric rotor experiment 
in a granular gas. The first state, with a Gaussian pdf 
for the angular velocity, is found to be well described 
by a first order Langevin equation and therefore analo- 
gous to a Brownian system of temperature T r . This is 
reflected in the observation that the FT is valid at large 
t for the injected work and for all t for the total entropy 
production. In the second state the pdf of uj is double- 
peaked due to convection rolls developing in the granular 
gas. Here symmetry relations are found for very large 
and very small r, of which the latter of leads to a good 
prediction of the ratio of the height of the two peaks. 
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[1] J. Kurchan, J. Stat. Mech., P07005 (2007) 

[2] D. J. Evans, E. G. D. Cohen, and G. P. Morriss, Phys. 

Rev. Lett. 71, 2401 (1993). 
[3] G. Gallavotti and E. G. D. Cohen, Phys. Rev. Lett., 74 

2694 (1995); G. Gallavotti and E. G. D. Cohen, J. Stat. 

Phys. 80, 931 (1995). 
[4] D. J. Evans and D. J. Searles, Adv. Phys. 51, 

(2002). 

[5] S. Ciliberto, S. Joubaud and A. Petrosyan, J. 

Mech., P12003 (2010). 
[6] K. Feitosa, N. Menon, Phys. Rev. Lett. 92, 164301 

(2004). 

[7] P. Eshuis, K. van der Weele, D. Lohse and D. van der 
Meer, Phys. Rev. Lett. 104, 248001 (2010) 

[8] J. Talbot, J. R. D. Wildman and P. Viot, Phys. Rev. Lett. 
107, 138001 (2011); J. Talbot, A. Burdeau and P. Viot, 



1529 



Stat. 



5 



J. Stat. Mech., P03009 (2011); J. Talbot, A. Burdeau and 
P. Viot, Phys. Rev. E 82, 011135 (2010). 
[9] R. Balzan et al, Phys. Rev. E 83, 031310 (2011). 
[10] G. Costantini, U. Marini Bettolo Marconi and A. Puglisi, 
Phys. Rev. E 75, 061124 (2007); G. Costantini, A. Puglisi 
and U. Marini Bettolo Marconi, J. Stat. Mech., P07004 
(2009). 

[11] B. Cleuren and R. Eichhorn, J. Stat. Mech., P10011 
(2008). 

[12] B. Cleuren and C. van den Broeck, Europhys. Lett. 77, 

50003 (2007). 
[13] U. Seifert, Phys. Rev. Lett. 95, 040602 (2005). 
[14] M. Esposito and C. Van den Broeck, Phys. Rev. Lett. 

104, 090601 (2010). 
[15] D. Lacoste and K. Mallick, Phys. Rev. E 80, 021923 



(2009). 

[16] This is the steady state FT. A variant is the detailed FT 
which holds for arbitrary time spans r. 

[17] Theoretically, the ratchet temperature is related to the 
granular temperature by T T — T g (l + e)/2. Both are or- 
ders of magnitude larger than ksT . 

[18] It should be stressed that here heat and entropy are not 
the thermal quantities but are defined by analogy to those 
for the stochastic motion of the rotor, i.e., with the tem- 
perature T r . 

analogous with the corresponding 
[19] Without the asymmetric coating, both directions are 
(within experimental and statistical errors) found to have 
the same probability. 



